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Memorie deiia 




Abstract. Using a chemical evolution model we investigate the intriguing suggestion that 
there are populations of stars in some globular clusters (e.g. NGC 2808, to Centauri) with 
enhanced levels of helium (Y ~ 0.28 to 0.40) compared to the majority of the population 
that presumably have a primordial helium abundance. We assume that a previous generation 
of massive low-metallicity Asymptotic Giant Branch (AGB) stars has polluted the cluster 
gas via a slow stellar wind. We use two independent sets of AGB yields computed from 
detailed models to follow the evolution of helium, carbon, nitrogen and oxygen in the cluster 
gas using a Salpeter initial mass function (IMF) and a number of top-heavy IMFs. In no 
case were we able to fit the observational constraints, Y > 0.30 and C+N+O ~ constant. 
Depending on the shape of the IMF and the yields, we either obtained Y > 0.30 and large 
increases in C+N+O or Y < 0.30 and C+N+O ~ constant. These results suggest that either 
AGB stars alone are not responsible for the large helium enrichment or that any dredge-up 
from this generation of stars was less than predicted by standard models. 
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AGB and post-AGB 



1. Introduction 

Understanding the history and evolution of 
galactic globular cluster (GC) stars poses one 
of the greatest challenges to astrophysics. The 
star-to-star abundance variations of the light 
elements C, N, O, Na, Mg and Al ob served 
in ev e ry well studie d cluste r to date ( iKraftl 
1 1994 iGratton et all 12004 and references 
therein) are not found in field stars of the 
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same metallicity ( Gratto rTet all 1200 0). Hence 
these abundance anomalies are somehow 
the result of the cluster environment. The 
variations of the elements follow a common 
pattern from cluster to cluster: C-N, O-Na and 
Mg-Al are all negatively correlated ilShetrone}. 
19961 iKrafteTatl fT997t ICannon et all Il998t 
Gratto netall I200U ICoh en & Melendez, 
20051 ICohenetall 120051) . The abundances 
of iron-peak, s and r-process elements do 
not show the same star-to-star scatter as the 
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light elements nor do these elements vary 
with the lig ht elem ents l|Gratton et ail 12004 
Uames et all l2004t lYongetall 120051) . The 
key points are that O has been destroyed in 
some stars by up to one order of magnitude, 
the C+N+O and Mg+Al abundances remain 
almost constant regardless of the absolute 
spread and there is no evidence for large-scale 
variation of neutron- capture elements, except 
in the case of to Cen ( Smi th et aTl l2~000). 

The self-pollution scenario, first proposed 
by ICottrell & Da Costal (I1981I) . is the most 
promising to explain these abundance trends. 
This is because the star-to-star abundance vari- 
ations in C, N, O and Na are observed in 
stars at or near the main-sequence turn-off 
in addition to stars along the giant branch 
llGratton et all200lHRamirez & Coherl 1 2003 1 : 
James et alll2004t ICohen & Melendell2005t 
Cohe n et all 120051) . The other scenario, deep 
mixing, is still required to operate in low-mass 
giants to conve rt some C to N a fter the lu- 
minosity bump (Charbonnel, 1994). Owing to 
the constant [Fe/H] of stars in a given GC it 
has been assumed that the source of the pollu- 
tion was intermediate-mass Asymptotic Giant 
Branch (AGB) stars with initial masses be- 
tween ~3 to 8M0 rather than supernovae. The 
hot bottom burning experienced by these ob- 
jects provides an environment (at least quali- 
tatively) to produce helium, c onvert C and O 
to N Ne to Na and Mg to Al ( lLattanzio et all 

2004) . The mass lost via the slow winds of 
AGB stars could, in principle, have been re- 
tained by the cluster from wh ich ne w stars may 
have been born dThoul et all 12002). 

There is an increasing amount of ev- 
idence for helium enrichment from hori- 
zontal branch (HB) morphology and main- 
seque n ce colour-magnitud e diagra ms (Norris, 
2004; ID'Antona & Caloit Eooi iLee et all 
2005; ICaloi & D'Antonal l2005t iPiotto et all 

2005 ) where the helium has been proposed to 
have come from a previous generation of low- 
metallicity intermediate-mass AGB stars. For 
example, the unusual HB morphology of NGC 
2808, which exhibits an extended blue tail 
and a gap separatin g the red and blue clumps 
(Bedi net all |2000), can be best explained if 
the blue stars have an enhanced amount of 



helium (up to Y ~ 0.32) compared to those 
in the red HB clum p wi t h a primordial Y ~ 

0. 24 (ID'Antona & Caloil 120041) . Further evi- 
dence for an enrichment in helium in NGC 
2808 comes fro m the p eculiar main sequence 
( ID'Antona et all I2005I) . where the bluer stars 
are inferred to have Y m 0.40. Observations 
by IPiotto et all ll2005l) showed that the blue 
main-sequence of to Centauri is more metal- 
rich than the red sequence, contrary t o what 
is exp ected from stellar evolution, and iNorrisI 
(2004) showed that isochrones with Y - 0.40 
fit the bluest stars. These intriguing pieces of 
observational evidence have motivated us to 
study the AGB self-pollution scenario from a 
global perspective. 

Here we use the Fenneretal. (2004) GC 
chemical evolution model to follow the evo- 
lution of helium in the intracluster gas. We 
explore AGB model uncertainties by using 
two independent sets of AGB yields, including 
those used in the previous study which were 
tailor made for NGC 6752 with a metallic - 
ity [Fe/H] « -1.4. This metallicity is slightly 
more metal-rich than the average metallicity of 
[Fe/H] ~ -1.6 for NGC 6752, to Centauri, M3 
an d Ml 3 with d ata taken from the catalogue 
of lHarrisI (Il996l) . The cluster NGC 2808 is 

more metal-rich with [Fe/H] 1.15 ( lHarrisL 

1996) but we feel these AGB yields are suit- 
able for this study because the yields would not 
change significantly in this [Fe/H] range. We 
also follow the evolution of C, N and O since 
they impose important empirical constraints, 

1. e. C+N+O « constant that must be met by 
the model. 

2. Helium Production in Asymptotic 
Giant Branch Stars 

The helium yields f rom the AGB stars used by 
iFenner et all ( 1200 4) are shown in table^as the 
total mass of 4 He expelled into the intraclus- 
ter medium by each mo del. We hereafter re- 
fer to the yields used by IFenner et all (2004) 
as our yiel ds. We also show the Z = 0.001 
yields from IVentura et all ll2002l) for compar- 
ison and no te that they agree to within ~ 30%. 
Our yields ( IFenner et all |20 04|) are sv s temat- 
ically larger even though Ven tura et al.1 ( 2002) 
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Table 1. The mass of helium (in Mq) expelled into the intracluster medium 



Model type 


1.25 


2.5 3.5 


5.0 5.5 6.5 


Fenner et al. (2004) 


0.160 


0.490 0.680 


1.53 - 1.92 


Ventura et al. C20021 




0.675 


1.14 1.30 



use a different convective model and mass-loss 
rate and observe shallow dredge-up in their 
computations. The relatively good agreement 
comes about because most of the helium mixed 
to the surface in m > 3Mq models comes 
from the second dredge-up, not the thermally- 
pulsing AGB phase. Also, helium yields are 
more robust to reaction rate uncertainties than 
other species (e.g. 23 Na) owing to the fact that 
the net result of hydrogen burning is helium fu- 
sion, regardless of the internal rates of the var- 
ious cycles (CNO, NeNa and MgAl). 

3. The initial mass function 

The globular cluster chemical evolution model 
was described in detail in Fenner et al. (2004). 
We made two main changes for this study: 
we u sed an ind ependent set of AGB yields 
from IVentura et all f2002), and we changed 
the shape of the in itial mass function; see 
iKarakas et all (120061) for more details. There 
are other AGB yields besides those we con- 
sider here although most are from synthetic 
computations (e.g. Izzard et al. (2004)) and for 
more metal-rich populations. Herwie (2004) 
has low-metallicity AGB yields computed 
from detailed models but the results of study 
would not change significantly if we were to 
have used these yields instead. This is because 
his results are somewhere between the Ventura 
et al. and the Monash models. That is they have 
deep dredge-up but only a few (less than 20) 
thermal pulses, wher eas the Mona s h mod els of 
~ 5Mq presented in Fenner et al. ( 2004J) have 
~ 100 thermal pulses plus deep dredge-up. 

In this section we focus on one crucial el- 
ement of the model - the initial mass function 
(IMF). Most discussions in the literature point 



towar d a universal mass function (|kroupa, 
120011) although the re are hints for va riations 
at low metallicity. iD'Antona & Caloil ( 120041) 
comment that to produce the amount of he- 
lium required to form the number of blue HB 
stars requires a factor of ~ 10 more 4 to 7Mq 
stars than produced by a standard Salpeter 
IMF. Stars between 1 to 3M0 are either not 
produced or ejected from the cluster owing to 
a lack of evidence for enhanced levels of s- 
process elements or carbon. However, there is 
little physical motivation why such an unusual 
IMF would result from the early evolution of 
GCs when there is no observational evidence 
that such strange IMFs formed in other systems 
of comparable mass such as dwarf spheroidal 
galaxies. Although the abundance patterns ob- 
served in these systems are remarkably dif- 
ferent from those in (most) globular clusters 
dPritzl et alll200l there is still no evidence for 
a top heavy IMF. 

With these considerations in mind, we 
compute simulations with a standard Salpeter- 
like IMF with a slope of 1.31, a flat Salpeter 
with slope of 0.30 and a top-heavy IMF that 
increases the number of intermediate-mass 
AGBs by a factor of 10, see fig-H We hereafter 
refer to the top-heavy IMF as "IMS-enhanced" 
(intermediate-mass star enhanced). 

4. Results 

In fig. 12 we show the temporal evolution 
of the C+N+O isotopes (on a log scale) 
as a function of the helium mass fraction 
assuming the flat Salpeter IMF. Using our 
yields we see a substantial increase in he- 
lium with a maximum abundance of Y » 
0.33, similar to (but smaller than) the value 
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Fig. 1. Possible initial mass functions for the first generation of globular cluster stars. We consider three 
possible choices for the IMF including a Salpeter with slope s = 1.31, a flat Salpeter with s = 0.30 and an 
IMF the IMS-enhanced IMF places about 10 times more mass in 3.5 to 6.5M© stars. 




0.40 
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Fig. 2. The temporal evolution of the C+N+O abundance ([CNO/Fe]) in the cluster gas as a function of 
Y assuming a flat S alpeter IMF with slope = 0.30. In the (left) panel we s how results using yields from 
lFenneretal1<2004) and in the (right) panel yields from Ventur aet alJl2002h . 
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required by isochrones llD' Antona & CaloiL 
12004 ID' Antona et ail 120051) to match the 
bluest HB and main-sequence stars of NGC 
2808. The maximum helium abundance is 
reached in under 200 Myr reflecting the dom- 
inant contribution of intermediate-mass AGB 
stars with lifetimes < 120 Myr. The increase 
in helium is also accompanied by a ~ 1 dex 
increase in C+N+O although there is some de- 
pletion of O by ~ 0.3 dex. 

The simulation using the Ventura et al. 
yields has a moderate increase in the total 
C+N+O of ~ 0.5 dex but this is probably 
within observational errors. The maximum Y 
in this case does not exceed 0.27, well below 
the inferred value of the bluest main sequence 
and HB stars. However, the substantial O de- 
pletion of ~ 0.5 dex in this case is similar 
to (but still smaller than) the maximum dis- 
persion observed in GC stars. If we compare 
to observations, the most "pollut ed" stars in 
M13 have [O/Fe] * -0.8 ( Kraft etail fl99l 
whereas "normal" stars have [O/Fe] +0.4, in- 
dicating significant O destruction of more than 
one order of magnitude. The simulation with 
the flat Salpeter produces many more stars with 
m > IMq compared to a standard Salpeter 
(see fig.Q resulting in more carbon along with 
helium injected into the intracluster medium. 
Indeed, this IMF produces results quite dif- 
ferent to the IMS-enhanced IMF which gives 
higher weight to the most massive AGB stars 
which have more efficient HBB and a lower 
weight to stars with 1 < ot(Mq) < 3. 
Results using the standard Sal peter and the 
IMS-e nhanced are presented in Karak aset alJ 
(2006) and are also not consistent with the ob- 
servational constraints we have considered in 
this study. 

5. Discussion & Conclusions 

Our investigation into the chemical evolution 
of helium in globular clusters has shown the 
difficulty the AGB self-pollution scenario suf- 
fers when trying to explain the large helium 
enrichment (Y > 0.30) hypothesized to fit 
the horizontal branch morphology of clusters 
like NGC 2808. We utilize a chemical evolu- 
tion model and two independent sets of AGB 



yields to follow the evolution of the intraclus- 
ter medium for a typical cluster with [Fe/H] 
= -1.4. We have tested three different IMFs 
for the first generation of GC stars including 
a standard Salpeter, a flat Salpeter and one 
IMF that boosts the number of intermediate- 
mass AGB stars by a factor of 10. The flat 
Salpeter produces maximum helium mass frac- 
tions in-between those found when using the 
standard Salpeter (with much lower Y values) 
and the IMS-enhanced I MF with slig htly larger 

Y abundances dKarakas et all |2006). The be- 
haviour of the C, N and O elements is however 
quite different, owing to the overall increased 
number of 1 to 6.5Mq stars which produce C 
as well as helium. 

Simulations assuming a flat Salpeter with 
a slope of 0.30 show helium mass fractions 
as large as Y ~ 0.33 but only with enor- 
mous increases in the total C+N+O content 
of the cluster gas, in violation of observations. 
The Ventura et al. yields predict a maximum 

Y ~ 0.27 with the total C+N+O abundance in- 
creasing by ~ 0.5 dex. Assuming that such an 
IMF is realistic we still have a problem fitting 
all of the observational constraints. Indeed, the 
use of such an IMF does not help the diffi- 
culties faced by the self-pollution scenario in 
matching the constraints that we have consid- 
ered in this study i.e. Y > 0.30 and C+N+O 
« constant. We not e that the previous study by 
iFenner et alJ ( l20 04) suffered similar difficulties 
when trying to match the observed abundances 
of the Na, Mg and Al isotopes in the cluster 
NGC 6752. 

There are many model uncertainties and 
in particular the extent of the dredge-up is 
far from known and shallower dredge-up, as 
observed in the Ventura et al. models, would 
help keep C+N+O constant. The nucleosyn- 
thesis associated with hot bottom burning is 
also dependent on the convective model (or 
the mixing-length parameter a for MLT mod- 
els) and varying this will i mpact the yields 
llVentura & D' Antona! |2005). However, more 
efficient convection leads to larger luminosities 
which would likely drive mass loss, leading to 
shorter AGB lifetimes and possibly smaller he- 
lium yields. 
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Thus far we have assumed that the large 
helium abundances inferred from theoretical 
isochrones are accurate. These determinations 
are model dependent and suffer from many of 
the same uncertainties that afflict AGB models 
(e.g. convection) and hence it is hard to gauge 
just how reliable these helium determinations 
are. Even if they are overestimates with the 
largest Y closer to 0.30 instead of 0.40, we still 
have a problem fitting the observational con- 
straints with the current set of AGB models. 
Perhaps given the difficulties associated with 
the self-pollution scenario we need to look to 
other solutions including pollution from out- 
side the clu ster, as discussed in th e case of co 
Centauri bv lBekki & Norrisl ( 120061) . 
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